Coherence freeze in an optical lattice investigated via pump-probe spectroscopy 
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Motivated by our observation of fast echo decay and a surprising coherence freeze, we have de- 
veloped a pump-probe spectroscopy technique for vibrational states of ultracold *^Rb atoms in an 
optical lattice to gain information about the memory dynamics of the system. We use pump-probe 
spectroscopy to monitor the time-dependent changes of frequencies experienced by atoms and to 
characterize the probability distribution of these frequency trajectories. We show that the inferred 
distribution, unlike a naive microscopic model of the lattice, correctly predicts the main features of 
the observed echo decay. 
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Characterizing decoherence mechanisms is a crucial 
task for experiments aiming to control quantum systems, 
e.g., for quantum information processing (QIP). In this 
work, we demonstrate how two-dimensional (2D) pump- 
probe spectroscopy may be extended to provide impor- 
tant information on these mechanisms. As a model sys- 
tem, we study quantum vibrational states of ultracold 
atoms in an optical lattice. In addition to being a leading 
candidate system for QIP [1] , optical lattices are proving 
a versatile testing ground for the development of quan- 
tum measurement and control techniques [2, 3] and a 
powerful tool for quantum simulations, e.g. the study of 
Anderson localization and the Hubbard model [4]. 

In our experiment, we study the vibrational coher- 
ence of ®^Rb atoms trapped in a shallow one-dimensional 
standing wave. Through our 2D pump-probe technique, 
we obtain detailed microscopic information on the fre- 
quency drift experienced by atoms in the lattice, en- 
abling us to predict the evolution of coherence. Since 
the pioneering development of the technique in NMR[5], 
2D spectroscopy has been widely used to obtain high- 
resolution spectra and gain information about relax- 
ations, couplings, and many-body interactions, in realms 
ranging from NMR [6] to molecular spectroscopy [7-10] 
to semiconductor quantum wells [11, 12]. Here, we show 
that similar powerful techniques can be applied to the 
quantized center-of-mass motion of trapped atoms, and 
more generally, offer a new tool for the characterization 
of systems in QIP and quantum control. 

We have previously measured the evolution of coher- 
ence between the lowest two vibrational states of poten- 
tial wells [14]. The dephasing time is about 0.3ms (T^). 
This dephasing is partly due to an inhomogeneous distri- 
bution of lattice depths as a result of the transverse Gaus- 
sian profile of the laser beams. To measure the homoge- 
neous decoherence time {T2), we perform pulse echoes, 
measuring the echo amplitude as a function of time [14]. 
Figure 1 shows two typical measurements of echo ampli- 
tude carried out on different dates under slightly different 
conditions such as different average lattice depths and 



different dephasing times. The echo amplitude initially 
decays with a time constant of about 0.7ms, which is 
much faster than the photon scattering time (~ 60ms) in 
the lattice. It then exhibits a l?7is-long coherence freeze 
followed by a final decay. Absent real decoherence on the 
short time scale of 1ms, only loss of frequency memory 
would inhibit the appearance of echoes. This loss comes 
about when atoms experience time-varying frequencies. 
We use 2D pump-probe spectroscopy to monitor this fre- 
quency drift. Our 2D pump-probe spectroscopy is essen- 
tially a version of spectral hole-burning for vibrational 
states. By monitoring the changes in the hole spectrum 
as a function of time we gain information on the atoms' 
frequency drift. Information obtained from our 2D spec- 
tra enables us to characterize the temporal decay of fre- 
quency memory and through our simulations we find that 
"coherence freeze" is related to the shape of this mem- 
ory loss function. Similar plateaus in echo decay and a 
two-stage decay of echo amplitude have been observed 
in a Cooper-pair box [16], for a single electron spin in 
a quantum dot [17] and for electron spins in a semicon- 
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Figure 1. (Color online) Two typical measurements of echo 
amplitude vs. time. The echo pulse and the observed echo 
envelope are centered at times tp and 2tp, respectively. Af- 
ter an initial decay, echo amplitude stays constant for about 
1ms forming a plateau, before decaying to zero. The average 
lattice depths are 20Er (circles) and ISEr (squares). 



ductor [18]. Those plateaus or two-stage decays have 
been either explained through a priori models or simply 
described phenomenologically. Here, we are introducing 
an experimental technique to directly probe the origin of 
plateaus. 

The periodic potential in our experiment is formed by 
interfering two laser beams blue-detuned by 25GHz from 
the D2 transition line, F = 3^i^' = 4(A = TSOnm), 
thus trapping atoms in the regions of low intensity, which 
minimizes the photon scattering rate and the transverse 
forces. The two laser beams intersect with parallel linear 
polarizations at an angle oi 9 — (49.0 ± 0.2)°, resulting 
in a spacing oi L ^ (0.930 ± 0.004)/im between the wells. 
Due to gravity, the full effective potential also possesses 
a "tilt" of 2.86i?fl per lattice site, where Er — g^^^ is 
the effective lattice recoil energy. The photon scatter- 
ing time in our experiment is « 60ms and the Landau- 
Zcnncr tunneling times for transitions from the lowest 
two levels are greater than 160tos. Atoms are loaded to 
the lattice during a molasses cooling stage and prepared 
in the ground vibrational state by adiabatic filtering [13]. 
Due to the short coherence length of atoms in optical mo- 
lasses (60nm at lOnK), there is no coherence between the 
wells. We measure populations of atoms in the ground 
vibrational, the first excited, and the (lossy) higher ex- 
cited states Pi, P2, and Pl, respectively, by fluorescence 
imaging of the atomic cloud after adiabatic filtering [13]. 

The pump and probe pulses are sinusoidal phase mod- 
ulations of one of the laser beams forming the lattice. 
The modulation is of the form (j){t) = A{t)[l — cos{u}mt)], 
where A(t) is a square envelope function with amplitude 
2it/72 and ujm is a variable frequency. The duration of 
each pulse is 8 cycles, i.e., T = 8{2tt/uj„i). This phase 
modulation shakes the lattice back and forth periodically, 
coupling vibrational states of opposite parity. To first- 
order in modulation amplitude, the phase-modulating 
part of the Hamiltonian has the same form as the electric 
dipole Hamiltonian. 

The inhomogeneous spectrum of vibrational excita- 
tions is measured in an average lattice depth of 2AEb, 
by applying probe pulses at different frequencies and 
measuring state populations. Figure 2(a) shows state 
populations Pi, P2, and Pl (black circles) as a func- 
tion of probe frequency. We then measure the pump- 
probe spectrum for a fixed delay. A pump pulse with 
a specific frequency is applied, exciting atoms in wells 
whose vibrational transition frequency matches the fre- 
quency of the pump pulse, therefore burning a hole in 
the spectrum of the ground state population. After a 
delay, probe pulses at different frequencies are applied, 
coupling the ground and excited states of atoms whose 
frequencies match those of the probe pulses. The red 
squares in Fig. 2(a) show populations Hi, H2, and Hi, 
as a function of probe frequency for a delay of 2ms and 
at a pump frequency of 6A5kHz. The central part of 
Hi shows the hole burnt into the spectrum of the ground 
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Figure 2. (Color online) Experiment: (a) black circles: probe- 
alone populations, Pi,P2, and Pl- Red squares: pump- probe 
populations, ni,Il2,nL for a pump at 6A5kHz and a delay 
of 2ms. The circled region shows the hole-burning signal in 
Pl and Hi. (b) The difference spectrum APi = Hi — Pi. (c) 
Growth in the r.m.s. width of a Gaussian fit to APi as a 
function of delay. The bar centered at (600 ± 25) Hz is the 
inhomogeneous width determined from the Pi spectrum. The 
green line centered at (417 ±5)^2 indicates the instrumental 
width expected due to the spectra of our pulses. 



state population. To characterize the hole, we plot the 
difference between the pump-probe and the probe-alone 
spectra, APi = Hi — Pi, in Fig. 2(b). We monitor the 
frequency drift of atoms by changing the delay between 
the pump and probe pulses. Figure 2(c) shows the r.m.s. 
width of APi as a function of delay. The r.m.s. width 
increases with increasing delay until it approaches the 
inhomogeneous width of the lattice. For pump-probe de- 
lays shorter than 2ms, the coherence present between 
the lowest two states results in Ramsey fringes making 
it impractical to extract useful spectra from the mea- 
surements. Repeating the above procedure for different 
pump frequencies, we obtain a 2D spectrum of APi for 
a fixed delay as a function of luq, the pump frequency at 
time io, and a;„, the probe frequency at i„. We inter- 
pret the measured 2D spectrum as the conditional prob- 
ability, P(a;„|a;o), of an atom having frequency a;„ at t„ 
given that at time to its frequency was ujq. The joint 
probability is P(w„,a;o) = P(a;„|a;o)P(a;o), where P(wo) 
is the probe-alone spectrum. Figures 3(a-b) show the ex- 
perimental joint probabilities P(a;2,a;o) and P(a;5,a;o), 
for delays of 2ms and 5ms, respectively. By measuring 
joint probabilities for three different delays, P(oj2,ujo), 
P{uj3,ujq) , and P(a;5,a;o), we construct a joint distribu- 
tion of three frequencies at three different times in order 
to characterize how frequencies change as a function of 
time (frequency trajectories) in the first 5ms. Details of 
our analysis are provided in the next two paragraphs. 

We assume a skew-normal distribution S9i{u;) [19], 
which is a generalization of the normal distribution. 
The multivariate skew-normal distribution is defined 
as 55v;(a;) = 2^(a;;fi)$(a-^a;), where ip is the three- 
dimensional normal density function, and ^ is the cu- 
mulative distribution function, u: = (tZ) — /x), where 
(D = (a;o,a;2,W5) denotes frequencies at three different 



times, /x = (/zq, /i2, A^s) and a = {ao, 012, as) are the three 
frequency means and shape parameters, respectively. 0, 
is the covariance matrix whose diagonal terms are vari- 
ances {erf}, and whose off-diagonal terms are covariances 
{atajPij}, where {pij} are the mutual correlation coeffi- 
cients. We estimate the best values for these parameters 
by using a genetic algorithm, as described below. 

The bivariate marginals of 55v;(a;), denoted 
as M(Lij„,ciJm), are calculated by integrating 
S9l with respect to the third frequency, e.g., 
Af(w2,'^o) — J duj5 S9i{uj). These marginals are 
then convolved with the power spectra of the pump 
and probe pulses to obtain the three bivariate con- 
volved marginals g{ujn,'-^m) for comparison to the three 
measured spectra, P{iOn,ujm)- Similarly, the univariate 
marginals of S9l{uj) are calculated and convolved with 
the power spectra of the probe pulses for comparison 
to the three measured probe-alone spectra, P(a;„). 
Then, the sum of the squared residuals for all data 
points is calculated. The convolved marginal 5 (0^5, 0^2) 
is compared to the experimental P{uj3,ujq) with the 
assumption that the frequency changes experienced by 
atoms in the (2 -^ 5)ms interval are the same as in the 
(0 -> 3)tos interval. From repeated runs of the genetic 
algorithm function in MatLab, we find 12 sets of best fit 
parameters, and take their average. Figures 3(c-f) show 
the two convolved marginals g{uj2,uJo), g{u!^,uJo) and 
the corresponding bare marginals M (0^2,^0), M{uj5,ujo) 
thus obtained. The distribution parameters are means 
{Mo,M2,M5} = {5.92(4), 5.92(3), 5.96(6)}/cifz; r.m.s. 



widths {CTo,o-2,a-5} = {0.77(5), 0.82(4), 0.83(7)}fcifz; 
mutual correlation coefficients {po2,P25,Po5} = 
{0.88(2), 0.82(3), 0.80(4)}; and shape parameters 
{ao,a2,a5} = {2.6(5), 3.2(4), 4.0(2.0)}. The numbers 
in parentheses are the standard deviations calculated 
from the 12 optimization runs. The means and r.m.s. 
widths stay approximately constant. Furthermore, the 
mutual correlations P25 and po5 are essentially equal 
and within two standard deviations of po2- The increase 
in the anti-diagonal width of M(a;5,a;o) compared to 
M{uj2,u}o) is a signature of the slight decrease in the 
mutual correlation coefficient pos compared to po2- 

The bare joint distribution s^K (t*^) is our estimate of 
the probability of frequency trajectories. Since s^ is 
a function of frequencies at three different times, fre- 
quency trajectories are constructed by interpolating be- 
tween frequencies at these times for all possible combi- 
nations of wq, UJ2, and ^5. Figures 4(a-c) show three 
groups of smoothed [20] trajectories that are constrained 
to originate from the frequency range Awq and to pass 
through a window of width Aa;2 a,t t = 2ms. Drawn 
for each group are small samples of trajectories selected 
according to their probability. The sideways curves 
are the final probabilities at t = 5ms calculated from 
Pf{uj^) — J^^ d(jjQ J^^ d(jj2S9i{u}). The most probable 
trajectories (Fig. 4(a)) stay within one standard devia- 
tion (width of the univariate skew-normal marginal at 
t = 0) from the center of distribution. The next most 
probable trajectories (Fig. 4(b)) are the ones where fre- 
quency initially increases in the first 2ms followed by 
trajectories where frequency decreases in the first 2ms 
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Figure 3. Contour plots: measured joint probabili- 
ties (a) P{u}2,i^o), (b) P(cj5,cJo); convolved marginals (c) 
(/(tJ2, <jJo), (d) g{iJ5,iJo) obtained by convolving bare marginals 
M{liJ„,ijJo) with the power spectra of the pump and probe 
pulses; bare marginals (e) M{uj2,ojo), (f) M{uj5,ujo)- 
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Figure 4. (Color online) (a,b,c) (cow from top left): Sam- 
ples of smoothed frequency trajectories selected according to 
their probabilities. Trajectories are constrained to originate 
from the frequency range Acjq = 0.5kHz and to pass through 
a window of width Alj^ — 0.52kHz at t — 2ms. Aujq is 
centered at 6.5kHz and Aa;2 is centered at: (a) 6.5kHz, (b) 
7.07kHz, and (c) 5.93kHz. The sideways curves show the 
final probabiUties Pf{uj5). (d) Projection of 55V;(a;) onto the 
plane defined by frequency differences (tJ2 —1^0) and (0)5 —^02), 
showing a negative correlation between the two time intervals. 



(Fig. 4(c)). Figure 4(d) is a projection of S9l{u;) onto 
the plane defined by frequency differences, {0J2 — wq) and 
(w5 — ^2)- Each quadrant corresponds to either an in- 
crease or decrease in the two time intervals, (0 -> 2)ms 
and (2 -> 5)ms. The negative correlation between the two 
time intervals indicates that a rising frequency in the first 
time interval is usually followed by a falling frequency in 
the second interval and vice versa. This is in sharp con- 
trast to the positive correlation we would obtain by con- 
sidering the ballistic expansion of atoms through a spatial 
Gaussian distribution of well depths. In that case, there 
are no trajectories where an initial frequency decrease 
is followed by an increase. Although ballistic expansion 
would give rise to similar time scales for echo amplitude 
decay it does not predict the formation of plateaus. This 
shows the power of 2D spectroscopy in characterizing all 
frequency trajectories without reference to any a priori 
model. 

The echo amplitude can be predicted by calculating the 
net accumulated phase, (/5(t ; u}), along each trajectory 
before and after the echo pulse, and by weighting each 
trajectory according to its probability S9l{u>), 



e(r) 



dui S9i{uj) 



^-iip{T;u>) 
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where ^(r ; O) = J^""'*^ dt"oj{t")-J^' dt'uj{t'). As Fig. 5 
shows, the decay time constant of echo amplitude pre- 
dicted by the bare probability distribution sc^l (w) is sim- 
ilar to the measured value. Depending on the exact 
method of interpolation (smooth or linear), the details 
of the echo decay vary but we always see a two-stage 
decay. We hypothesize that the non-decaying values of 
mutual correlation coefficients P25 and po5 lead to the 
formation of the plateau. By repeating simulations with 
reduced values of P25 and pos: we find that for values of 
(/O25 — Po5 > 0.4) the plateau disappears. A more accu- 
rate calculation of echo decay would require knowledge of 
the joint probability distribution of multiple frequencies 
at multiple times on a finer time scale. The complexity 
of the task of obtaining data and extracting a joint dis- 
tribution grows rapidly as the time scale is made finer. 
The extracted SfK. i<^) is our best estimate for a trivariate 
joint distribution. An exact reconstruction of the dis- 
tribution requires knowledge of higher-order conditional 
probabilities such as P(w5|a;2,a;o)- In principle, it would 
be possible to extract this conditional probability from 
a two-pump, one-probe measurement. Simulations, how- 
ever, show that this signal might be very difficult to ex- 
tract in practice [21]. 

In conclusion, we have developed a two-dimensional 
pump-probe spectroscopy technique for vibrational states 
of ultracold atoms in optical lattices to gain information 
about the memory dynamics of the system without any 
need for a priori physical models. Our method is a di- 
rect frequency-domain characterization of the system for 




Figure 5. (Color online) A selection of echo curves, ^(t), gen- 
erated by different s^H. (t*') resulting from repeated optimiza- 
tion trials. The inset shows the echo decay curve generated 
by S!^l (w) calculated from the average values of parameters. 



different delays between the pump and probe pulses. By 
measuring populations of the quantum vibrational states, 
we infer information about the underlying correlations, 
which in turn enables us to predict how coherence evolves 
in the system. By directly extracting a best estimate 
trivariate joint distribution s:H_ (w) from the measured 2D 
spectra at three different delays, and constructing trajec- 
tories accordingly, we have calculated echo amplitude as 
a function of time. The inferred distribution correctly 
predicts the main features of the observed echo decay 
and it suggests that persistence of frequency-frequency 
correlations leads to the formation of plateaus. The sim- 
ilarity between the predicted echo decay curve and the 
directly measured ones demonstrates the power of 2D 
spectroscopy in providing information on memory dy- 
namics, information which could not always have been 
foreseen. Our 2D pump-probe method is a general tech- 
nique that should be broadly applicable in a variety of 
quantum information candidate systems for characteri- 
zation of decoherence mechanisms. 
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